Transcriptome and single-cell analysis reveal disulfidptosis-related modification patterns of tumor microenvironment and prognosis in osteosarcoma.

Osteosarcoma (OS) is the most common malignant bone tumor with high pathological heterogeneity. Our study aimed to investigate disulfidptosis-related modification patterns in OS and their relationship with survival outcomes in patients with OS. We analyzed the single-cell-level expression profiles of disulfidptosis-related genes (DSRGs) in both OS microenvironment and OS subclusters, and HMGB1 was found to be crucial for intercellular regulation of OS disulfidptosis. Next, we explored the molecular clusters of OS based on DSRGs and related immune cell infiltration using transcriptome data. Subsequently, the hub genes of disulfidptosis in OS were screened by applying multiple machine models. In vitro and patient experiments validated our results. Three main disulfidptosis-related molecular clusters were defined in OS, and immune infiltration analysis suggested high immune heterogeneity between distinct clusters. The in vitro experiment confirmed decreased cell viability of OS after ACTB silencing and higher expression of ACTB in patients with lower immune scores. Our study systematically revealed the underlying relationship between disulfidptosis and OS at the single-cell level, identified disulfidptosis-related subtypes, and revealed the potential role of ACTB expression in OS disulfidptosis.

SET domain containing 2 promotes megakaryocyte polyploidization and platelet generation through methylation of α-tubulin.

BACKGROUND: Megakaryocytes (MKs) are polyploid cells responsible for producing ∼1011 platelets daily in humans. Unraveling the mechanisms regulating megakaryopoiesis holds the promise for the production of clinical-grade platelets from stem cells, overcoming significant current limitations in platelet transfusion medicine. Previous work identified that loss of the epigenetic regulator SET domain containing 2 (SETD2) was associated with an increased platelet count in mice. However, the role of SETD2 in megakaryopoiesis remains unknown. OBJECTIVES: Here, we examined how SETD2 regulated MK development and platelet production using complementary murine and human systems. METHODS: We manipulated the expression of SETD2 in multiple in vitro and ex vivo models to assess the ploidy of MKs and the function of platelets. RESULTS: The genetic ablation of Setd2 increased the number of high-ploidy bone marrow MKs. Peripheral platelet counts in Setd2 knockout mice were significantly increased ∼2-fold, and platelets exhibited normal size, morphology, and function. By knocking down and overexpressing SETD2 in ex vivo human cell systems, we demonstrated that SETD2 negatively regulated MK polyploidization by controlling methylation of α-tubulin, microtubule polymerization, and MK nuclear division. Small-molecule inactivation of SETD2 significantly increased the production of high-ploidy MKs and platelets from human-induced pluripotent stem cells and cord blood CD34+ cells. CONCLUSION: These findings identify a previously unrecognized role for SETD2 in regulating megakaryopoiesis and highlight the potential of targeting SETD2 to increase platelet production from human cells for transfusion practices.

Differentiation route determines the functional outputs of adult megakaryopoiesis.

Emerging evidence has revealed a direct differentiation route from hematopoietic stem cells to megakaryocytes (direct route), in addition to the classical differentiation route through a series of restricted hematopoietic progenitors (stepwise route). This raises the question of the importance of two alternative routes for megakaryopoiesis. Here, we developed fate-mapping systems to distinguish the two routes, comparing their quantitative and functional outputs. We found that megakaryocytes were produced through the two routes with comparable kinetics and quantity under homeostasis. Single-cell RNA sequencing of the fate-mapped megakaryocytes revealed that the direct and stepwise routes contributed to the niche-supporting and immune megakaryocytes, respectively, but contributed to the platelet-producing megakaryocytes together. Megakaryocytes derived from the two routes displayed different activities and were differentially regulated by chemotherapy and inflammation. Our work links differentiation route to the heterogeneity of megakaryocytes. Alternative differentiation routes result in variable combinations of functionally distinct megakaryocyte subpopulations poised for different physiological demands.

Metabolic Landscape of Osteosarcoma: Reprogramming of Lactic Acid Metabolism and Metabolic Communication.

BACKGROUND: Lactic acid, previously regarded only as an endpoint of glycolysis, has emerged as a major regulator of tumor invasion, growth, and the tumor microenvironment. In this study, we aimed to explore the reprogramming of lactic acid metabolism relevant to osteosarcoma (OS) microenvironment by decoding the underlying lactic acid metabolic landscape of OS cells and intercellular signaling alterations. METHODS: The landscape of OS metabolism was evaluated using single-cell gene expression data, lactic acid metabolism clustering, and screening of the hub genes in lactic acid metabolism of OS samples using transcriptome data. The role of the hub gene NADH:Ubiquinone Oxidoreductase Complex Assembly Factor 6 (NDUFAF6) was validated with in vitro studies and patient experiments. RESULTS: Single-cell RNA sequencing data validated a lactic acid metabolismhigh subcluster in OS. Further investigation of intercellular communications revealed a unique metabolic communication pattern between the lactic acid metabolismhigh subcluster and other subclusters. Next, two lactic acid metabolic reprogramming phenotypes were defined, and six lactic acid metabolism-related genes (LRGs), including the biomarker NDUFAF6, were screened in OS. In vitro studies and patient experiments confirmed that NDUFAF6 is a crucial lactic acid metabolism-associated gene in OS. CONCLUSIONS: The patterns of lactic acid metabolism in OS suggested metabolic reprogramming phenotypes relevant to the tumor microenvironment (TME) and identified NDUFAF6 as an LRG prognostic biomarker.

Single-cell sequencing and transcriptome analyses in the construction of a liquid-liquid phase separation-associated gene model for rheumatoid arthritis.

Background: Rheumatoid arthritis (RA) is a disabling autoimmune disease that affects multiple joints. Accumulating evidence suggests that imbalances in liquid-liquid phase separation (LLPS) can lead to altered spatiotemporal coordination of biomolecular condensates, which play important roles in carcinogenesis and inflammatory diseases. However, the role of LLPS in the development and progression of RA remains unclear. Methods: We screened RA and normal samples from GSE12021, GSE55235, and GSE55457 transcriptome datasets and GSE129087 and GSE109449 single-cell sequencing datasets from Gene Expression Omnibus database to investigate the pathogenesis of LLPS-related hub genes at the transcriptome and single cell sequencing levels. Machine learning algorithms and weighted gene co-expression network analysis were applied to screen hub genes, and hub genes were validated using correlation studies. Results: Differential analysis showed that 36 LLPS-related genes were significantly differentially expressed in RA, further random forest and support vector machine identified four and six LLPS-related genes, respectively, and weighted gene co-expression network analysis identified 396 modular genes. Hybridization of the three sets revealed two hub genes, MYC and MAP1LC3B, with AUCs of 0.907 and 0.911, respectively. Further ROC analysis of the hub genes in the GSE55457 dataset showed that the AUCs of MYC and MAP1LC3B were 0.815 and 0.785, respectively. qRT-PCR showed that the expression of MYC and MAP1LC3B in RA synovial tissues was significantly lower than that in the normal control synovial tissues. Correlation analysis between hub genes and the immune microenvironment and single-cell sequencing analysis revealed that both MYC and MAP1LC3B were significantly correlated with the degree of infiltration of various innate and acquired immune cells. Conclusion: Our study reveals a possible mechanism for LLPS in RA pathogenesis and suggests that MYC and MAP1LC3B may be potential novel molecular markers for RA with immunological significance.

DPP3 expression promotes cell proliferation and migration in vitro and tumour growth in vivo, which is associated with poor prognosis of oesophageal carcinoma.

Dipeptidyl peptidase III (DPP3), a zinc­dependent metallopeptidase, is upregulated in a variety of malignancies. However, little is known about its roles in the pathogenesis of these malignancies. The present study was designed to investigate the roles of DPP3 in the pathogenesis and progression of oesophageal cancer (EC). The expression level of DPP3 in EC tissues and adjacent normal tissues was detected in 93 cases of tissue biopsies collected from patients diagnosed with oesophageal carcinoma by immunohistochemistry. The effect of DPP3 expression on cell proliferation, migration or apoptosis was determined in DPP3­depleted EC cells created by infection with lentivirus containing short hairpin RNA specific to the human DPP3 mRNA sequence, followed by detection at the cellular level using a Celigo cell count assay, flow cytometry, wound­healing assay and Transwell assay as well as chip screening with a Human Apoptosis Antibody Array kit, which enables the quantitative detection of 43 apoptosis­related genes. A xenograft model was applied to detect the tumour growth and invasion of DPP3­depleted cancer cells in nude mice. The results revealed that DPP3 expression was elevated in EC tissues compared with adjacent non­tumour tissues, and high DPP3 expression was significantly associated with poor prognosis. DPP3 depletion resulted in reduced cell proliferation and migration and enhanced cell cycle arrest and apoptosis of EC cells and led to the inhibition of tumour growth and invasion in a xenograft model. In addition, DPP3 depletion was associated with the upregulation of the proapoptotic proteins SMAC and p53 and the downregulation of the antiapoptotic proteins clAP­2, IGFBP­2 and TRAILR­4. Finally, DPP3 may promote cell proliferation, migration and survival of EC cells in vitro and tumour growth and invasion of oesophageal carcinoma in vivo, and thus may serve as a molecular target for tumour therapy.

NBAS, a gene involved in cytotoxic degranulation, is recurrently mutated in pediatric hemophagocytic lymphohistiocytosis.

Hemophagocytic lymphohistiocytosis (HLH), particularly primary HLH (pHLH), is a rare, life-threatening disease. Germline genetic deficiency of 12 known HLH genes impairs cytotoxic degranulation in natural killer (NK) cells or cytotoxic T lymphocytes (CTLs) and contributes to pHLH development. However, no pathogenic mutations in these HLH genes are found in nearly 10% of HLH patients, despite a strong suspicion of pHLH, suggesting that the underlying genetic basis of HLH is still unclear. To discover novel susceptibility genes, we first selected 13 children with ppHLH (presumed primary HLH patients in the absence of detectable known HLH gene variants) and their parents for initial screening. Whole-genome sequencing (WGS) in one trio and whole-exome sequencing (WES) in twelve trios revealed that two ppHLH patients carried biallelic NBAS variants, a gene that is involved in Golgi-to-endoplasmic reticulum (ER) retrograde transport upstream of the degranulation pathway. Additionally, two candidate genes, RAB9B and KLC3, showed a direct relationship with known HLH genes in protein-protein interaction (PPI) network analysis. We analyzed NBAS, RAB9B, KLC3 and known HLH genes in an independent validation cohort of 224 pediatric HLH patients. Only biallelic NBAS variants were identified in three patients who harbored no pathogenic variants in any of the known HLH genes. Functionally, impaired NK-cell cytotoxicity and degranulation were revealed in both NBAS biallelic variant patients and in an NBAS-deficient NK-cell line. Knockdown of NBAS in an NK-cell line (IMC-1) using short hairpin RNA (shRNA) resulted in loss of lytic granule polarization and a decreased number of cytotoxic vesicles near the Golgi apparatus. According to our findings, NBAS is the second most frequently mutated gene (2.11%) in our HLH cohort after PRF1. NBAS deficiency may contribute to the development of HLH via a dysregulated lytic vesicle transport pathway.

Single-cell Transcriptomic Analysis Reveals the Cellular Heterogeneity of Mesenchymal Stem Cells.

Ex vivo-expanded mesenchymal stem cells (MSCs) have been demonstrated to be a heterogeneous mixture of cells exhibiting varying proliferative, multipotential, and immunomodulatory capacities. However, the exact characteristics of MSCs remain largely unknown. By single-cell RNA sequencing of 61,296 MSCs derived from bone marrow and Wharton's jelly, we revealed five distinct subpopulations. The developmental trajectory of these five MSC subpopulations was mapped, revealing a differentiation path from stem-like active proliferative cells (APCs) to multipotent progenitor cells, followed by branching into two paths: 1) unipotent preadipocytes or 2) bipotent prechondro-osteoblasts that were subsequently differentiated into unipotent prechondrocytes. The stem-like APCs, expressing the perivascular mesodermal progenitor markers CSPG4/MCAM/NES, uniquely exhibited strong proliferation and stemness signatures. Remarkably, the prechondrocyte subpopulation specifically expressed immunomodulatory genes and was able to suppress activated CD3+ T cell proliferation in vitro, supporting the role of this population in immunoregulation. In summary, our analysis mapped the heterogeneous subpopulations of MSCs and identified two subpopulations with potential functions in self-renewal and immunoregulation. Our findings advance the definition of MSCs by identifying the specific functions of their heterogeneous cellular composition, allowing for more specific and effective MSC application through the purification of their functional subpopulations.

Prognostic implications of N6-methyladenosine RNA regulators in breast cancer.

The significance of N6-methyladenosine (m6A) RNA modifications in the progression of breast cancer (BC) has been recognised. However, their potential role and mechanism of action in the tumour microenvironment (TME) and immune response has not been demonstrated. Thus, the role of m6A regulators and their downstream target gene components in BC remain to be explored. In this study, we used a series of bioinformatics methods and experiments to conduct exploratory research on the possible role of m6A regulators in BC. First, two regulatory modes of immune activation and inactivation were determined by tumour classification. The TME, immune cell infiltration, and gene set variation analysis results confirmed the reliability of this pattern. The prognostic model of the m6A regulator was established by the least absolute shrinkage and selection operator and univariate and multivariate Cox analyses, with the two regulators most closely related to survival verified by real-time quantitative reverse transcription polymerase chain reaction. Next, the prognostic m6A regulator identified in the model was crossed with the differential copy number of variant genes in invasive BC (IBC), and it was determined that YTHDF1 was a hub regulator. Subsequently, single-cell analysis revealed the expression patterns of m6A regulators in different IBC cell populations and found that YTHDF1 had significantly higher expression in immune-related IBC cells. Therefore, we selected the intersection of the BC differential expression gene set and the differential expression gene set of a cell line with knocked-down YTHDF1 in literature to identify downstream target genes of YTHDF1, in which we found IFI6, EIR, and SPTBN1. A polymerase chain reaction was conducted to verify the results. Finally, we confirmed the role of YTHDF1 as a potential prognostic biomarker through pan-cancer analysis. Furthermore, our findings revealed that YTHDF1 can serve as a new molecular marker for BC immunotherapy.

Characteristics of the Intestinal Microorganisms in Middle-Aged and Elderly Patients: Effects of Smoking.

Introduction: Smoking affects the occurrence and development of many diseases. We attempt to study the structure of intestinal flora in the middle-aged and elderly population as well as how smoking affects the intestinal flora. Methods: We collected population information, biochemical indicators, and patient feces from 188 middle-aged and elderly male patients, and their feces were tested for the 16S rRNA gene of intestinal flora. Results: We performed a cluster analysis on the intestinal structure of the included population and found that there was a significant difference in the number of smokers between each group (p = 0.011). Subsequently, the microbiological diversity analysis of current smokers and nonsmokers was carried out. The results indicated that there was a significant difference in species composition between the two groups (p = 0.029). Through the analysis on LEfSe differential bacteria, it was found that in current smoking patients, the abundances of the genus Bifidobacterium and the genus Coprobacillus were less, while the abundances of the genera Shigella, Paraprevotella, Burkholderia, Sutterella, Megamonas, and p-75-a5 under the family level of Erysipelotrichaceae were slightly high. We analyzed the correlation between the abundances of these eight different bacteria and clinical indicators. The results revealed the following: the abundance of the genus Bifidobacterium was negatively correlated with fasting blood glucose (r = -0.198, p = 0.006) and positively correlated with uric acid (r = 0.207, p = 0.004) and total bilirubin (r = 0.175, p = 0.017); Shigella bacteria were positively correlated with fasting blood glucose (r = 0.160, p = 0.028) and uric acid (r = 0.153, p = 0.036) levels; the genus Paraprevotella and BMI (r = -0.172, p = 0.018) are negatively correlated; the abundance of the genus Burkholderia was positively correlated with γ-glutamyltransferase (r = 0.146, p = 0.045) levels; Sutterella was correlated with fasting blood glucose (r = 0.143, p = 0.05) and creatinine level (r = -0.16, p = 0.027), which was positively correlated with fasting blood glucose and negatively correlated with creatinine. Conclusions: In middle-aged and elderly patients with cardiovascular disease, smoking can reduce the abundance of Bifidobacterium, while the abundances of some negative bacteria such as Burkholderia, Sutterella, and Megamonas increase.

A prospective study revealing the role of an immune-related eRNA, WAKMAR2, in breast cancer.

Enhancer RNAs (eRNAs) are a subclass of non-coding RNAs that are generated during the transcription of enhancer regions and play an important role in tumourigenesis. In this study, we focused on the crucial eRNAs that participate in immune responses in invasive breast cancer (IBC). We first used The Cancer Genome Atlas and Human enhancer RNA Atlas to screen for tissue-specific eRNAs and their target genes. Through Pearson correlation analysis with immune genes, the eRNA WAKMAR2 was identified as a key candidate involved in IBC. Our further research suggested that WAKMAR2 is crucial in regulating the tumour microenvironment and may function by regulating immune-related genes, including IL27RA, RAC2, FABP7, IGLV1-51, IGHA1, and IGHD. Quantitative reverse transcription-polymerase chain reaction was used to detect the expression of WAKMAR2 in IBC and normal tissues, and the effect of WAKMAR2 on the regulation of downstream genes in MB-231 and MCF7 cells was studied in vitro. WAKMAR2 was found to be highly involved in tumour immunity and was downregulated in IBC tissues. Furthermore, the expression of WAKMAR2 and its target genes was observed at the pan-cancer level. This study provides evidence to suggest new potential targets for the treatment of breast cancer.

Revealing the Immune Infiltration Landscape and Identifying Diagnostic Biomarkers for Lumbar Disc Herniation.

Intervertebral disc (IVD) degeneration and its inflammatory microenvironment ultimately led to discogenic pain, which is thought to originate in the nucleus pulposus (NP). In this study, key genes involved in NP tissue immune infiltration in lumbar disc herniation (LDH) were identified by bioinformatic analysis. Gene expression profiles were downloaded from the Gene Expression Omnibus (GEO) database. The CIBERSORT algorithm was used to analyze the immune infiltration into NP tissue between the LDH and control groups. Hub genes were identified by the WGCNA R package in Bioconductor and single-cell sequencing data was analyzed using R packages. Gene expression levels were evaluated by quantitative real-time polymerase chain reaction. The immune infiltration profiles varied significantly between the LDH and control groups. Compared with control tissue, LDH tissue contained a higher proportion of regulatory T cells and macrophages, which are associated with the macrophage polarization process. The most significant module contained three hub genes and four subclusters of NP cells. Functional analysis of these genes was performed, the hub gene expression pattern was confirmed by PCR, and clinical features of the patients were investigated. Finally, we identified TGF-ß and MAPK signaling pathways as crucial in this process and these pathways may provide diagnostic markers for LDH. We hypothesize that the hub genes expressed in the specific NP subclusters, along with the infiltrating macrophages play important roles in the pathogenesis of IVD degeneration and ultimately, disc herniation.

Pan-cancer analysis reveals tumor-associated macrophage communication in the tumor microenvironment.

BACKGROUND: Tumor-associated macrophages (TAMs) are abundant in the tumor microenvironment (TME). However, their contribution to the immunosuppressive status of the TME remains unclear. METHODS: We integrated single-cell sequencing and transcriptome data from different tumor types to uncover the molecular features of TAMs. In vitro experiments and prospective clinical tests confirmed the results of these analysis. RESULTS: We first detected intra- and inter-tumoral heterogeneities between TAM subpopulations and their functions, with CD86+ TAMs playing a crucial role in tumor progression. Next, we focused on the ligand-receptor interactions between TAMs and tumor cells in different TME phenotypes and discovered that aberrant expressions of six hub genes, including FLI1, are involved in this process. A TAM-tumor cell co-culture experiment proved that FLI1 was involved in tumor cell invasion, and FLI1 also showed a unique pattern in patients. Finally, TAMs were discovered to communicate with immune and stromal cells. CONCLUSION: We determined the role of TAMs in the TME by focusing on their communication pattern with other TME components. Additionally, the screening of hub genes revealed potential therapeutic targets.

Identification of Immune-Related Therapeutically Relevant Biomarkers in Breast Cancer and Breast Cancer Stem Cells by Transcriptome-Wide Analysis: A Clinical Prospective Study.

Cancer stem cells (CSCs) represent a subset of tumor cells that are responsible for recurrence and metastasis of tumors. These cells are resistant to radiotherapy and chemotherapy. Immunotherapeutic strategies that target CSCs specifically have provided initial results; however, the mechanism of action of these strategies is unclear. The data were requested from The Cancer Genome Atlas and Genotype-Tissue Expression, followed with the survival analysis and weighted gene co-expression network analysis to detect survival and stemness related genes. Patients were divided into three groups based on their immune status by applying single sample GSEA (ssGSEA) with proven dependability by ESTIMATE analysis. The filtered key genes were analyzed using oncomine, GEPIA, HPA, qRT-PCR, and functional analysis. Patients in a group with a higher stemness and a lower immune infiltration showed a worse overall survival probability, stemness and immune infiltration characteristics of breast cancer progressed in a non-linear fashion. Thirteen key genes related to stemness and immunity were identified and the functional analysis indicated their crucial roles in cell proliferation and immune escape strategies. The qRT-PCR results showed that the expression of PIMREG and MTFR2 differed in different stages of patients. Our study revealed a promising potential for CSC-target immunotherapy in the early stage of cancer and a probable value for PIMREG and MTFR2 as biomarkers and targets for immunotherapy.

Effects of the cooking modes on commonly used pesticides residue in vegetables and their chronic dietary exposure risk in South China.

Effects of cooking modes on the real intake and chronic exposure risk of pesticide residues in vegetables are usually neglected and largely unknown. Four modes of daily meal preparation; chafing dish, soup, salad and stir-frying were studied in this work to clarify their impact on the residual pesticides in foods. A detection method for 14 types of pesticide residues in different cuisines was developed. In this work, chronic exposure risks of four pesticides were analysed by probabilistic assessment based on data from public health and a pesticide residues investigation conducted. The results showed that chafing dish and soup methods greatly lowered the types, contents and exposure risks from residue pesticides. Salad preserved almost all the pesticide residues, and the risks were also relatively high in the stir-frying method. In chafing dish and soup, pesticide residues were dispersed in the media and posed quite low threats to humans. Considering the age, infants and children were at a higher risk of exposure than other populations. Reassuringly, all of the risks were at acceptable levels. This study clarified how the cooking modes affect chronic exposure risks to pesticide residues in the vegetables. The outcomes also show the effects of cooking method on healthy daily diets.

Construction and analysis of a spinal cord injury competitive endogenous RNA network based on the expression data of long noncoding, micro­ and messenger RNAs.

Spinal cord injury (SCI) results from trauma and predominantly affects the young male population. SCI imposes major and permanent life changes, and is associated with high future mortality and disability rates. Long non-coding RNAs (lncRNAs) have recently been demonstrated to serve critical roles in a broad range of biological processes and to be expressed in various diseases, including in SCI. However, the precise mechanisms underlying the roles of lncRNAs in SCI pathogenesis remain unexplored. In the present study, the aim was to identify critical differentially expressed lncRNAs in SCI based on the competing endogenous RNA (ceRNA) hypothesis by mining data from the Gene Expression Omnibus database of the National Center for Biotechnology Information and to unveil the functions of these lncRNAs. Different approaches and tools were employed to establish a network consisting of 13 lncRNA, 93 messenger RNA and 9 microRNA nodes, with a total of 202 edges. Three node lncRNAs were identified based on the degree distribution of the nodes, and their corresponding subnetworks were subsequently constructed. Based on these subnetworks, the biological pathways and interactions of these 3 lncRNAs were detailed using FunRich software (version 3.0). The primary results of the 3 lncRNA enrichment analyses were that they were associated with autophagy, extracellular communication and transcription factor networks, respectively. The phosphoinositide 3­kinase/protein kinase B/mammalian target of rapamycin signaling pathway of XR_350851 was the classic autophagy pathway, indicating that XR_350851 may regulate autophagy in SCI. The possible role of XR_350851 in SCI revealed in the current study based on the regulatory mechanism of ceRNAs has uncovered a new repertoire of molecular factors with potential as novel biomarkers and therapeutic targets in SCI.

Application of network construction to estimate functional changes to insulin receptor substrates 1 and 2 in Huh7 cells following infection with the hepatitis C virus.

Hepatitis C virus (HCV) is closely associated with insulin resistance (IS), acting primarily by interfering with insulin signaling pathways, increasing cytokine-mediated (tumor necrosis factor α, interleukin 6) inflammatory responses and enhancing oxidative stress. In the insulin signaling pathways, the insulin receptor substrate (IRS) is one of the key regulatory factors. The present study constructed gene regulatory sub­networks specific for IRS1 and IRS2 in Huh7 cells and HCV­infected Huh7 (HCV­Huh7) cells using linear programming and a decomposition algorithm, and investigated the possible mechanisms underlying the function of IRS1/2 in HCV­induced IS in Huh7 cells. All data were obtained from GSE20948 of the Gene Expression Omnibus database from the National Center for Biotechnology Information. Genes which were significantly differentially expressed between Huh7 and HCV­Huh7 cells were analyzed using the significance analysis of microarray algorithm. The top 50 genes, including IRS1/2, were used as target genes to determine the gene regulatory networks and next the sub­networks of IRS1 and IRS2 in HCV­Huh7 and Huh7 cells using Gene Regulatory Network Inference Tool, an algorithm based on linear programming and the decomposition process. The IRS1/2 sub­networks were divided into upstream/downstream groups and activation/suppression clusters, and were further analyzed using Molecule Annotation System 3.0 and Database for Annotation, Visualization, and Integrated Discovery software, two online platforms for enrichment and clustering analysis and visualization. The results indicated that in Huh7 cells, the downstream network of IRS2 is more complex than that of IRS1, indicating that the insulin metabolism in Huh7 cells may be primarily mediated by IRS2. In HCV­Huh7 cells, the downstream pathway of IRS2 is blocked, suggesting that this may be the underlying mechanism in HCV infection that leads to insulin resistance. The present findings add a further dimension to the understanding of the pathological mechanisms of HCV infection-associated insulin resistance, and provide novel concepts for insulin resistance and glucose metabolism research.

Computational networks of activating transcription factor 3 gene in Huh7 cell lines and hepatitis C virus-infected Huh7 cell lines.

Activating transcription factor 3 (ATF3) is an adaptive­response gene of the ATF family. ATF3 activity may be induced in response to a number of different stress-associated signals and ATF3 is involved in a variety of cellular processes. However, the functions of ATF3 and its molecular networks in human hepatoma cells lines and hepatitis C virus-infected Huh7 (HCV-Huh7) cells are not well understood. In the present study, ATF3 regulatory networks in Huh7 and HCV-Huh7 cell lines were established using the linear programming-based GRNinfer software and molecule annotation system 3.0 software. The gene expression omnibus dataset, GSE20948, was analyzed. The resulting network consisted of clusters located upstream and downstream of ATF3 in Huh7 and HCV-Huh7 cell lines. Using the annotation, visualization and integrated discovery (DAVID) software, 10 activation and 2 inhibition enriched functional annotation clusters were identified downstream of ATF3 in HCV-Huh7 cells. However, there were no enriched functional annotation clusters identified upstream of ATF3 in HCV-Huh7 cells. Furthermore, no clusters were identified downstream nor upstream of ATF3 in Huh7 cells. Gene ontology term and Kyoto encyclopedia of genes and genomes pathway analyses demonstrated that ATF3 may be involved in a number of biological processes, in particular, in metabolism regulation in HCV-Huh7 cells. It is hypothesized that the ATF3 pathway may be activated in Huh7 cells following HCV infection and that it is a potential 'hub' in the network of HCV-Huh7 cells.

hOGG1 Ser326Cys polymorphism and risk of hepatocellular carcinoma among East Asians: a meta-analysis.

BACKGROUND: The hOGG1 gene encodes a DNA glycosylase enzyme responsible for DNA repair. The Ser326Cys polymorphism in this gene may influence its repair ability and thus plays a role in carcinogenesis. Several case-control studies have been conducted on this polymorphism and its relationship with the risk of hepatocellular carcinoma (HCC) among East Asians. However, their results are inconsistent. METHODS: We performed a meta-analysis of published case-control studies assessing the association of the hOGG1 Ser326Cys polymorphism with HCC risk among East Asians. PubMed, EMBASE, SCI, BIOSIS, CNKI and WanFang databases were searched. A random-effect model was used to calculate odds ratios (ORs) and 95% confidence intervals (95% CIs). Analyses were conducted for additive, dominant and recessive genetic models. RESULTS: Eight studies were identified involving 2369 cases and 2442 controls assessing the association of the hOGG1 Ser326Cys polymorphism with HCC risk among East Asians. Applying a dominant genetic model, only in the Chinese population, the Cys allele was significantly associated with increased risk of HCC (OR 1.56, 95% CI 1.12-2.17). However, two studies influenced this finding according to sensitivity analysis. Furthermore, considerable heterogeneity and bias existed among Chinese studies. CONCLUSION: There is limited evidence to support that the hOGG1 Ser326Cys polymorphism is associated with HCC risk among East Asians. Well-designed and large-sized studies are required to determine this relationship.